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Abstract 

Aims. We try to understand the gas heating and cooling in the S 140 star forming region by spatially and spectrally resolving the dis¬ 
tribution of the main cooling lines with GREAT/SOFIA and combining our data with existing ground-based and Herschel observations 
that trace the energy input and the density and temperature structure of the source. 

Methods. We mapped the fine structure lines of [Oi] (63 pm) and [Cii] (158 pm) and the rotational transitions of CO 13-12 and 
16-15 with GREAT/SOFIA and analyzed the spatial and velocity structure to assign the emission to individual heating sources. We 
measure the optical depth of the [C ii] line and perform radiative transfer computations for all observed transitions. By comparing the 
line intensities with the far-infrared continuum we can assess the total cooling budget and measure the gas heating efficiency. 

Results. The main emission of fine structure lines in S 140 stems from a 8.3" region close to the infrared source IRS 2 that is not 
prominent at any other wavelength. It can be explained by a photon-dominated region (PDR) structure around the embedded cluster if 
we assume that the [O i] line intensity is reduced by a factor seven due to self-absorption. The external cloud interface forms a second 
PDR at an inclination of 80-85 degrees illuminated by an UV field of 60 times the standard interstellar radiation field. The main 
radiation source in the cloud, IRS I, is not prominent at all in the fine structure lines. We measure line-to-continuum cooling ratios 
below lO^'', i.e. values lower than in any other Galactic source, rather matching the far-IR line deficit seen in ULIRGs. In particular 
the low intensity of the [Cii] line can only be modelled by an extreme excitation gradient in the gas around IRS 1. We found no 
explanation why IRS I shows no associated fine-structure line peak, while IRS 2 does. 

Conclusions. The inner part of S 140 mimics the far-IR line deficit in ULIRGs thereby providing a template that may lead to a future 
model. 

Key words. ISM: individual (S 140) - ISM: structure - ISM: clouds - ISM: photon-dominated region (PDR) - ISM: lines and bands 
- ISM: abundances 


1. Introduction 

The cooling budget of the gas in galaxies is dominated by the 
emission in the fi ne structure lines of atomic oxygen and ionized 
carbon (see e.g. iBurton et al.l [l990l: iRollig et alJ l2007l) . Their 
emission stems to a large fraction from UV-illuminated surfaces 
of molecular clouds, so called photon-dominated regions ( PDRs 
iTielens & Hollenbachlll985l : ISternberg & Dalgarnoill995h . Due 

to the clumpy or fractal structure of molecular clouds, the PDR 
surfaces may represen t a large fraction of their total volume 
(lOssenkopf et alj2007l) . To quantify the importance of the PDRs 
in terms of the overall cooling balance of molecular clouds, we 
have to compare the fine structure line emission with the lines 
of CO (and its isotopologues) that measure the emission of the 
molecular material and to the continuum emission from interstel¬ 
lar dust. The dust emission measures the combination of the to¬ 
tal column density structure with the local heating, optically thin 
molecular lines provide us with a view on the spatial and veloc¬ 
ity distribution of the molecular gas, and the fine structure lines 
contain the information on the local energy input through UV ra¬ 
diation and the cooling balance in the gas. The German REceiver 
for Astronomy at Terahertz Erequencies (GREAT) onboard the 
Stratospheric Observatory for Infrared Astronomy (SOFIA) al¬ 
lows for the first time to obtain velocity resolved spectra of 


the 63 pm line of a tomic oxygen, [Oi], with a good sensitivity 
(iBiichel et alj|2015l) . We used the instrument to measure the line 
simultaneously with the 158 um line of ionized carbon, [C ii], in 
the S 140 star-forming regioiQ. 

S 140 is an Hii-region at the surface of the L 1204 molec¬ 
ular cloud, created by the illuminat ion from the BOV star 
HD 211880, at a distance of 764 pc (iHirota et al.ll200^ . The 


external surface of the m olecular cloud is exposed to a mod- 
erate UV field of 140yn (iPoelman & SDaanil200-5l) to 230;^fo 
(iTimmerm ann et all 1 19961) (xo = standard interstellar radiation 


field [S 


irame 


1978h . However, ongoing star-formation within 


the cloud alre ady formed a clus ter of deeply embedded radi¬ 
ation sources (lEvans et al.l Il989h . probably creating more in¬ 
ternal PDRs around the sources. Radio continuum observa¬ 
tions detected ultracompact H n region s at the strongest sources 
(e.g. iTofani et~an Il995t [Hoarel l2()06l) . The strongest submm 
source, IRS 1, harbors a cluster of massive young stellar ob¬ 
jects at a projected distance of 75" from the cloud inter- 


' Comparison to integrated line observations from the Herschel 
Space Observatory are provided in Appx. \M 

Herschel is an ESA space observatory with science instruments pro¬ 
vided by European-led Principal Investigator consortia and with impor¬ 
tant participation from NASA. 
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face (iMinchin et alJ Il993h . IWeigelt et al.l (l2002l) showed the 
presence of an outflow cavity around IRSl. The correspond¬ 
ing outflow walls are illuminate d by the cluster so that they 
form internal high-density PDRs. IPedes et al] (l2010ll suggested 
that the [C ii] emission from IRS 1 measured by HlFl onboard 
Herschel stems from these irradiated outflow walls, exposed 
to much higher radiation fields of more than lO^^o- Physical 
models of the struc t ure around IRS 1 have been proposed b y 
Harvey et al.l(ll978h.lGurtler et al. (11991 ). IMinchin et al.l(ll993h. 


^^derTM.et alJ (2000)7 iMueller et alJ0002h. Ide Wit et all 

( 20091) . iMaud et alJ ( 2013). The whole cloud is fragmented 
and clumpy. [Poelman & Snaansl (l2006l) determined densities of 
nHj ~ 4 X 10^ cm~^ for the clumps and nnj ~ 10^ cm“^ for the 
interclump medium. 

In Sect. 2, we present the new GREAT observations of the 
63 //m [Oi] and 158 jum [Cii] fine structure and the CO 16-15 
and 13-12 lines towards S 140. In Sect. 3, we discuss the mor¬ 
phology and velocity structure of the observed emission com¬ 
pared to other lines, in particular low-/ CO transitions, and lit¬ 
erature data. In Sect. 4, we derive gas parameters across the map 
and for selected positions to understand the nature of the emit¬ 
ting gas. Sect. 5 discusses possible reasons for the association of 
the fine structure emission peak with an otherwise inconspicuous 
source and the overall extremely low fine-structure line emis¬ 
sion. 


2. GREAT observations 

We used the German REceiver for Astronom y at Terahertz 
Erequencies (GREAT0, iHevminck et al.l l2Q12h onboard the 
Stratospheric Obs ervatory for Infrared Astronomy (SOFIA, 
I Young et al.ll2012h to take on-the-fly maps of the central region 
ofS 140 in [Oi],63//m, [Cii],CO 16-15 and CO 13-12. 

The [Oi] observations were performed in the nights from 
May 16 and 17, 2014. The [Cii] data were taken on May 16, 
17, and 20, 2014. CO 13-12 was observed in two outer strips on 
May 20, 2014 and CO 16-15 was measured in the inner strip on 
January 13, 2015. 

We mapped an area of 144 x 96" oriented at an angle of 37° 
east of north, perpendicular to the outer interface of the PDR. 
Due to instrument scheduling constraints, the [Oi] and CO 16-15 
observations were restricted to an inner strip of only 40" width 
covering the central cluster, while the observation of the CO 
13-12 line avoided the central 32" covering only the outer two 
strips. In [Cii] we have the complete map covering the whole 
area. All observations involving [Oi] were taken on a 3" grid 
to obtain full sampling, the outer strips were mapped with an 
8" sampling, guaranteeing a close-to full sampling for the tele¬ 
scope beam at the [C ii] frequency. The OFF position for the ob¬ 
servations was located south-west of H ii region and molecular 
cloud at R.A.= 22''18”37^0, Dec= 63°14'18".0 (J2000), six ar- 
cminutes away from IRS 1, our reference position in the maps at 
R.A.= 22^'19”18^21, Dec= 63°18'46".9. XFFT spectrometers 
with a resolution of 44 kHz were used as backends. 

Atmospheric calibration was d one through the standard 
GREAT pipeline dGuan et al.ll2012h . The pointing accuracy of 
the telescope and the alignment of the GREAT instrument rela¬ 
tive to the telescope optical axis give a combined pointing error 
of less than 2 arcseconds. The two GREAT channels operating 
in parallel were co-aligned to 1 arcsec accuracy. The forward 

^ GREAT is a development by the MPI fr Radioastronomie and 
the KOSMA / Universitat zu Koln, in cooperation with the MPI fr 
Sonnensystemforschung and the DLR Institut fr Planetenforschung 


efficiency is 0.97 for all observed frequencies. The other param¬ 
eters of the GREAT beam at the observed frequencies are given 
in Tab. [1] All spectra were scaled to main beam Rayleigh-Jeans 
temperatures. 

Eor the spectra of CO and [C ii] the default pipeline provided 
quite flat baselines. To be independent of the continuum level, 
a first order baseline, measured over 100 kms * outside of the 
15 kms^' of the line was subtracted. As the [Oi] maps suffered 
from some OTF-striping, these data were further processed by a 
secondary OEF subtraction, treating the first and last points of all 
OTF lines as reference allowing for the correction of gain drifts. 
This assumes that the map edges fall in regions free of emission, 
an assumption that is confirmed when inspecting the raw data. 
To improve the signal to noise ratio, the spectra were smoothed 
with a box-car Alter to a resolution of 0.3-0.4 km s '. Due to 
the different combinations of observations, the total noise in the 
sum spectra taken within one beam is variable over the map. The 
average rms noise levels of the individual spectra and the beam- 
averaged sum spectra across the map at the smoothed resolution 
are given in Tab. [1] 

Some of the [O i] spectra showed more irregular, wavy base¬ 
line structures not described by a simple linear relation. We tried 
to correct them by fitting various orders of polynomials to the 
baseline varying the order from two to ten. We did not And a sig¬ 
nificant global improvement of the spectra with increasing order, 
but an improvement of a few individual spectra. By using the 
criterion of the least negative values in the integrated intensity 
map, we selected the seventh order for all spectra but use the 
outcome of the variation of the baseline orders basically as an 
error estimate for the baseline uncertainty. Eor different baseline 
orders, the intensities of individual bins in the line changed by 
up to 0.4 K and the integral over the full velocity range of the 
line (-15 - 0 km s ') varied by up to 6 K km s '. This uncer¬ 
tainty has to be added to the noise and the calibration error of the 
instrument of approximately 20 % as the uncertainty of the data. 

3. Spatial structure of the emission 

Figure [T] shows the observed integrated intensity map of the 
[C It] line superimposed on a Spitzer IRAC image of the cloud 
composed of the 3.6, 5.6 and 8 /tm channels. The IRAC image 
emphasizes two structures in the cloud. The peak in all bands 
is given by the embedded source IRS 1, relatively deep in the 
cloud. The embedded cluster produces about 10^ solar lumi- 
nosities, Lp, heating the dust to temperatures of up to 1400 K 
(IKoumnia et al.ll2015h . In the south-west we see the cloud sur¬ 
face illuminated by HD 211880 as a red band-like structure due 
to the excitation of PAHs at this external PDR. On top of these 
dominant features, we see some smaller spots east of IRS 1 close 
to IRS 3, bright in 5.6 emission. The [Cii] emission traces 
the front of the PDR in the south-west through a secondary peak, 
but has a clear maximum 20" north of IRS 1, slightly north of 
IRS 2. 

To compare the fine structure lines with the molecular cloud 
material we use the maps of CO 1-0, 2-1, C'^O 1-0, and '^CO 
1-0 taken with the IRAM 30m telescope, presented and dis¬ 
cussed by iKoumpia et akl (l2015h . Figure compares the [Cii] 
peak intensities with the CO 2-1 peak intensity map. The spatial 
distribution of the CO 2-1 line shows a behavior similar to the 
mid-infrared continuum. CO 2-1 peaks about 20" south of the 
[Cii] peak in a region covering IRS 1 and IRS 3. The CO 2-1 
maximum falls close to IRS 3 and at the cloud surface the line 
is brighter than deep in the cloud. There is a clear layering with 
[Cii] peaking closer to the illuminating star and CO 2-1 deeper 
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Table 1. Parameters of the GREAT line observations. 


Species 

Transition 

Vo 

[MHz] 

E„p 

[K] 

«crit“ 

[cm“3] 

beam EWHM 
[arcsec] 

7mb 

Av 

[kms^'] 

(T- * 

C' /,rms 

[K] 

(T ^ 

c' map,rms 

[K] 

C+ 

[Cii] "P 3/2 Pi/2 

1900536.90 

91.21 

3 X 103 

14.1 

0.67 

0.3 

3.7 

0.8-1.2 

O 

[0I]3P, -3^2 

4744777.49 

227.71 

5x 10^ 

6.6 

0.67 

0.4 

3.2 

1.3-1.9 

CO 

13-12 

1496922.91 

503.10 

1 X 10« 

18.6 

0.70 

0.3 

1.9 

0.9-1.3 


16-15 

1841345.51 

751.70 

2x lO*- 

14.5 

0.69 

0.4 

2.9 

0.7-1.0 


Notes. Collision rates from IWiesenfeld & Goldsmitl] 120 14h ; ITaguet et aO dl992l) : lYang et al.l 120101 at 50K). For optically thick lines the 
effective critical densities can be considerably lower due to radiative trapping. Noise in individual spectra. Noise in Nyquist-sampled map 
with GREAT beam FWHM. 
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Figure 1. Integrated intensity map of [C ii] overlaid on a Spitzer 
IRAC false-color image of the same region. The IRAC map uses 
the logarithmic intensities of the 3.6, 5.6 and 8//m channels. 
The integrated intensity contours are drawn at 50, 70, 90, 120, 
150, 180 K km s“*. The external cloud interface is visible in the 
south-west. The arrow shows the direction towards the illumi¬ 
nating source HD 211880. The red marks indicate the position 
of the known infrared sources: south-west is IRS 1, the north¬ 
ern o ne is IRS 2, and the south-east we find IRS 3 (lEvans et al.l 
Il989ll. the yellow marks show position of known submm cores 
iMinchin et al.lll9^lMaud et al.ll^OOh . The IRAC peak falls at 
the location of IRS 1. 



in the m olecular cloud, as expected f rom standard PDR models 
(see e.g. [Hollenbach & Tielerislll999h . The location of the peak 
of [C It] emission at the cloud surface therefore represents an ex¬ 
pected morphology, but the relative shift of the emission peak to 
the north of IRS 1 is a striking unexpected feature in the [C ii] 
maps. 

Eigure [3] compares the [Cii] with the [Oi] 63 //m line inte¬ 
grated intensity. The [O i] map shows a good match of the emis¬ 
sion peak with that of [Cii]. The peak is strongly concentrated 
and dominates the overall map. [Oi] also peaks slightly north 
of IRS 2; it may be offset from the center of the [C ii] peak by 
about 2", but this is close to our pointing accuracy. The location 
of the [C It] and [O i] peak indicates an association with IRS 2, 
similar to an already know n small continuum feature visible at 
mid-infrared wavelengths (iHarvev et al.ll201^ . The [Oi] peak 
is slightly elongated towards the south-east in the direction to¬ 
wards IRS 3. IRS 1 is not prominent at all in [O i]. We also see a 
slightly increased [O i] intensity at the cloud interface in the re- 
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Figure 2. Peak intensity map of [C ii] overlaid on a IRAM 30m 
map of the CO 2-1 peak intensity. The [C ii] contours are drawn 
at 15, 20, 25, 30, 35, 40, 45 K. The black marks indicate the 
positions of IRS 1-3. 
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Figure 3. Integrated intensity map of the [O i] 63 line (col¬ 
ors) overlaid by contours of the [C ii] integrated intensity (levels 
as in Fig. [1]). The black marks indicate the positions of IRS 1- 
3. All coordinates are taken relative to the position of IRS 1 at 
R.A.= 22*■19“18^21, Dec= 63°18'46".9 (J2000) 
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Figure 4. Integrated intensity map of CO 16-15 with contours of 
CO 13-12drawn at levels of 10,20, 30,40, 50, 60,70 K kms '. 
Only the region around IRS 1 has been observed in CO 16-15; 
in CO 13-12 only two outer strips have been covered avoiding 
the center. The black marks indicate the positions of IRS 1-3. 

gion where [C ii] peaks. Unfortunately, the data quality does not 
allow to discriminate whether there is any relative offset between 
the two fine structure line peaks in that region. 

Figure |4] shows the integrated intensities of the CO 13-12 
and 16-15 lines. Unfortunately, only complementary regions 
were mapped in the two lines with narrow overlapping strips. 
Therefore, the CO 13-12 observations are hardly useful to un¬ 
derstand the nature of the peak seen in the fine-structure lines. 
We see, however, that both high-7 CO lines show a more ex¬ 
tended emission around the central sources compared to the fine- 
structure lines. The CO 16-15 peak falls between IRS 1 and 
IRS 2, i.e. we see a superposition of hot CO emission from both 
IRS 1 and from the fine-structure line peak close to IRS 2. The 
high-7 CO lines trace both components in a similar way; they 
show an intermediate behavior between the fine-structure lines 
and the low-7 CO. 

3.1. Line profiles 

In all the published maps of S 140 only the '^CO 1-0 line peak 
intensity also has a global maximum close to IRS 41 To under¬ 
stand the nature of the fine-structure line peak it is, therefore, 
worth to compare our observations with the '^CO 1-0 data more 
carefully, including the velocit y structure. *^CO 1-0 i s optically 
thin over most parts of the map. lKoumnia et £0(12015h measured 
optical depths of about 0.3 for most of the gas and a maximum 
value of 0.8 that is still marginally optically thin. Figure|5]com- 
pares the distribution of the peak intensity, the mean line veloc¬ 
ity, and the line width of the ^^CO 1-0 line with the cotTespond- 
ing parameters from the [C ii] line. 

The relative proximity of the ^^CO 1-0 peak and the [C ii] 
and [O i] peaks seem to be accidental as the overall morphology 

^ For simplicity we use the location of the fine-structure line peak 
at about (AR.A.,A Dec) = (0",20") and IRS 2 synonymously in the 
following as the observed spatial deviations fall within the pointing ac¬ 
curacy of about 2". 


of the line maps shows no other similarities. The *^CO 1-0 emis¬ 
sion is much more extended and has about the same level at the 
positions of IRS 2 and IRS 1. When inspecting the '^CO 1-0 in¬ 
tegrated intensity, we find the same behavior as seen in the other 
molecular line tracers with a maximum close to IRS 1. A second 
similarity in the peak intensities of 1-0 and [C ii] are the 
relatively low values south of IRS 1. The spot match es the south¬ 
ern pa rt of the molecular outflow discussed e.g. by iMaud et al.l 
(l2013h . The outflow is driven from IRS 1 pointing towards the 
north-west and south-east, leading to the velocity gradient seen 
e.g. in the '^CO 1-0 line position map around IRS 1. At the loca¬ 
tion of the southern outflow [C ii] is hardly present at all, but the 
'^CO 1-0 line is prominent, tracing the outflow through broad 
line wings (see Fig. |5};) that lead to a large integrated intensity 
there in spite of the lower peak intensity. 

The [Cii] line shows the largest line width north-west of 
IRS 1, in the direction of the northern outflow, but only a moder¬ 
ate line width at the location of the intensity peak close to IRS 2. 
We And a trend of anti-correlation between the [C ii] lines and 
the '^CO 1-0 lines. The '^CO 1-0 lines are broad in a ridge 
south-east and north-east of the cluster that is associated with 
low center velocities while [C ii] shows the broadest lines west 
of the cluster. Except for the southern outflow region and the 
north-eastern boundary of the map, the [C ii] profiles are always 
somewhat broader than the '^CO 1-0 profiles. We interpret this 
as being due to optical depth broadening (see below). 

When inspecting the line positions, the '^CO 1-0 profiles 
are dominated by the outflow pattern with slightly blue shifted 
velocities south-east of IRS 1 and red shifted velocities in the 
opposite directiorQ. Overall, the variations are relatively small. 
In contrast, [C ii] shows a clear offset in the line velocity for the 
bright emission component close to IRS 2. While in most of the 
map the [C ii] line has the same velocity of about -8 km s ' as 
the bulk of the '^CO gas, the peak source has a significantly 
lower velocity of about -6.5 km s '. This indicates that the 
bright emission source is only weakly associated with the rest 
of the gas and not seen in '^CO. 

Figure |6] shows the line profiles of [Cii], [Oi], CO 16-15, 
CO 2-1, C'^O 1-0 and '^CO 1-0 at the location of the fine- 
structure line peak close to IRS 2, at the position of the main 
heating source IRS 1, and at the external cloud interface. All 
lines trace the interface at a velocity of about -8 km s '. The 
[Oi] line is only marginally detected there, [Cii] and the CO 
lines are narrow with a width of 3-4- km s '. The velocity of 
this PDR surface is slightly different from the large scale veloc¬ 
ity held of the thin surrounding medium measured through the 
broad H i line in the 36' beam of the Leid en/Dwingeloo survey 
to peak at -7 km s ' (iKalberla et al.ll200^ . 

Towards IRS 2, [Cii] peaks at -6.5 kms“', offset in veloc¬ 
ity from the bulk of the gas, as seen already in Fig. |5] The [O i] 
line shows a dip at this velocity and two peaks at the velocities 
of -5.5 km s ' and -8 km s ' coinciding with the peak veloc¬ 
ities seen in CO 2-1 towards IRS 2 and towards the interface. 
However, as there is only weak [O i] emission from the inter¬ 
face, we think that this coincidence is accidental and we rather 
see in [Oi] the same material that we trace in [Cii], indicating 
that [O i] exhibits a self-absorption dip where [C ii] peaks instead 
of two velocity components. 

CO 2-1 traces mainly other gas. We And bro ad outflow 
wings from the molecular outflow (IMaud et'^l2013h . only seen 


'' Follow-up investigations should resolve whether the broad, blue- 
shifted '^CO emission north-east of the cluster traces a second, so far 
unknown molecular outflow. 
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Figures. Peak intensity, first moment (centroid velocity) and second moment (translated to FWHM) of the lines of '^CO 1-0 
(colors) and [C ii] contours. The black marks indicate the positions of IRS 1-3. The contour levels for the [C ii] peak intensity match 
those from Fig.|2] for the line position they are given at -9.5, -8.5, -7.5, -6.5 kms ', and for the line width at 3, 4, 5, 6 kms“'. 
Velocities are given in LSR and intensities in main beam temperatures. 
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Figure 6 . Line profiles of the [Cii], [Oi], CO 16-15, CO 2-1, C'^O 1-0, and '^CO 1-0 lines at the position of the [Oi] peak (0,20 
close to IRS 2), at IRS 1 (0,0), and at the interface (-40,-54). The lines are averaged in a Gaussian beam of 17.3" FWHM (telescope 
-I- convolution kernel) except for the '^CO and C'*0 lines that can be only treated at the lower resolution of 23". 


in the low-7 CO lines. The optically thin lines of CO 16-15 and 
^^CO 1-0 are narrow at the peak suggesting that even the [Cii] 
line is broadened by a significant optical depth. The compari¬ 
son of the fine structure lines and CO 16-15 at IRS 1 shows 
that [C ii] only traces the velocity component matching the in¬ 
terface velocity there, but no emission at -6 kms ', while CO 
16-15 and [Oi] show a stronger emission in the -6 kms"^ com¬ 
ponent and a weaker emission with the interface velocity there. 
This means that we find several gas components with different 
velocities and chemical properties towards IRS 1. 

To estimate the optical depth of the [C ii] line, based on the 
line profiles, we can use two approaches; derive the optical- 
depth line broadening by comparing the line width with an op¬ 
tically thin tracer, and to compare the line intensity with an op¬ 
tically thin tracer not affected by abundance or excitation un¬ 
certainties, namely the [*^Cii] transitions. In the first approach, 
we assume a Gaussian velocity dispersion and obtain the opti¬ 


cal depth at the line peak from the measured line broadening 
o~iinp. (1 -H 0.1 15 t) cTypi where f is the line -center optical depth 
(iPhillips et al.lll97^IOssenkopf et alJl2013h . Using the optically 
thin 1-0 to measure the true velocity dispersion, we can 

translate Fig. |5}; into a [C ii] optical depth map, except for the 
region south of IRS 1 where the '^CO 1-0 line is broader than 
[C ii] due to the molecular outflow. For the interface, we obtain 
a 1-0 line width of 3 km s“^ and a [C ii] line width of 
4 km s ', in the region between IRS 1 and IRS 2, the line widths 
are 5 and 6 kms *, respectively. This translates into [Cii] peak 
optical depths of three towards the interface and two towards the 
central cluster. 

Figure Q shows the corresponding [*^Cii] spectra measured 
in the two regions. They are compared with a scaled version of 
the [C ii] line, reduced in intensity by the typ ical abun¬ 

dance ratio o f 1/67 in the solar neighborhood (iLanger & Penzia^ 
Il990l Il993h and the relative weights of the hyperfine compo- 
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Figure?. Baseline of the [Cii] profiles showing the hyperfine 
components of the ['^C ii] transitions. The upper plots shows the 
sum spectrum from all spectra with integrated intensities above 
100 K km s ' from the IRS 1-3 region, the lower plot from the 
interface region. The blue lines represent the same spectra scaled 
by an abundance ratio of 1/67 and the corresponding line ratio 
and frequency of the three hyperfine transitions, i.e. they repre¬ 
sent the spectrum that we expect in case of an optically thin [C ii] 
line. 


nents (0.625 at 11 km s“', 0.25 at -65 km s“', 0.125 at 63 km s“* 
relative to [Cii], lOssenkonf et al.ll2013h . Towards the interface 
the two stronger hyperfine components are clearly detected. The 
peak intensity of the [*^Cii] is 3.1 times as high as expected 
from the [Cii] line when assuming optically thin [Cii] emis¬ 
sion. This can be translated into a line-center optical depth of 
about three. Towards the embedded cluster, the strongest com¬ 
ponent appears as a shoulder of the broader [C ii] line and the 
two weaker components are barely visible. The factor 2.2 differ¬ 
ence between observed intensity and scaled optically thin [Cii] 
emission is also in agreement with a line-center optical depth 
of two as measured from the line broadening. Both approaches, 
therefore provide matching results for the [C ii] optical depth to¬ 
wards the two peaks of [C ii] emission. 

To resolve the intrinsic velocity profile of the main emission 
source at IRS 2 we compare scaled versions of the different pro¬ 
files towards this source in Fig. [8] The scaling factors in the plot 
were adjusted to allow for a comparison of the red wing of the 
profiles at sufficiently high velocities to not be affected by self¬ 
absorption or saturation from high optical depths. For CO 2-1 



V [km/s] 


Figures. Line profiles of the [Cii], [Oi], CO 2-1, CO 16-15, 
and *^CO 1-0 lines at the position of the [Oi] peak (0,20). The 
lines were scaled to similar intensities at -4 km s“*. For CO 2-1, 
a constant emission of 20 K is subtracted to exclude the extended 
outflow emission. 


we also had to subtract the very extended outflow contribution 
(see Fig. |6ll. The C'^O 1-0 line is not shown here as it matches 
almost exactly the shape of the CO 16-15 line. A comparison 
of the blue wing of the profiles provides no information on the 
source at IRS 2 as this velocity interval is heavily affected by 
emission and (self-) absorption from the bulk of the gas in S 140 
at -8 kms“*. 

We find an almost perfect match of the red wing shape of the 
three (partially) optically thick lines of [Cii], [Oi], and CO 2-1. 
In contrast, the thin lines of '^CO 1-0 and CO 16-15 are clearly 
narrower, not allowing to fit the same wing. Optical depth broad¬ 
ening of the lines must play a significant role. The peak of the 
CO 16-15 line coincides with the absorption dip seen in [Oi] 
and CO 2-1. The narrow shape of the peak of the CO 2-1 line at 
-4.7 km s ' results from this absorption out of a broader emis¬ 
sion component. The shoulder seen in the '^CO 1-0 and CO 2-1 
lines at -8 kms"' probably traces emission from IRS 1 or the 
extended molecular cloud material. Overall, the lines are consis¬ 
tent with the velocity dispersion as measured by the CO 16-15 
and C'^O 1-0 lines of only 2.4 km s“' (FWHM), but an optical- 
depth broadening of the fine-structure lines. 

For the optically thick lines, the intrinsic line intensity from 
the emission peak therefore should be higher than our observed 
integrated intensity. For [C ii], the profiles are consistent with the 
optical depth of about two, leading to a 45 % reduction of the in¬ 
tegrated line intensity compared to optically thin emission. For 
[Oi], we see, however, a clear self-absorption signature due to 
gas in the foreground with lower excitation temperature. That 
makes an estimate of the intrinsic intensity from the emission 
peak difficult. If we assume that the red wing is not affected by 
self absorption, we can use the ratio of the integrated line inten¬ 
sity between [Cii] and [Oi] seen in Fig. [8] as an estimate for the 
missing flux in the [Oi] line. This indicates that the source intrin¬ 
sic line flux should be higher by about 50 % without foreground 
absorption. 

From the information given in the line profiles we can al¬ 
ready conclude that the [O i] and [C ii] emission does not seem 
to be related to the known outflow activity in S 140, probably 
rather tracing PDR interfaces. While [Cii] and [Oi] are both 
bright at IRS 2, only [C ii] is bright at the cloud surface and [O i] 
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Figure 9. CO column density (upper plot; colors), gas density 
(lower plot: colors), and kinetic temperature (both plots: con¬ 
tours) derived from the RADEX fit of the integrated intensities 
of the observed lines of CO 1-0, 2-1, 16-15 and the isotopo- 
logues '^CO 1-0 and C'^O 1-0. 


is barely dete cted there. As the [Oi] line has a critical density of 
5x10^ cm“^ (I.Taauet et alJl992i see Tablll), 100 times the critical 
density of the [C ii] line, this probably reflects the lower density 
of the interface region. In contrast, at the fine structure line peak 
position, the gas must be dense and optically thick for [O i] so 
that the same material produces an almost Gaussian line in [C ii] 
and a double-peak-shaped line in [Oi]. The optically very thick 
CO 2-1 line only reveals a prominent red shifted wing towards 
IRS 2. This would be consistent with an expansion scenario of 
the source. 


4. Properties of the gas 

4.1. Gas parameters across the GREAT map 

As a first step to understand the nature of the fine structure 
lines we need to know the properties of the gas across the map. 
This can be done by analyzing the molecular lines of CO and 
its isotopologues based on our high-7 CO observations and the 
complementary low-7 CO, '^CO, and C'^O IRAM data from 
iKoumnia et al.l (l2015h . When knowing the column density, the 
volume density and kinetic temperature of the gas we can com¬ 
pute column density maps of [C ii] and [O i] from the measured 
fine-structure line intensities. 


We performed a fit of the CO (and isotopologues) line in- 
tensities using the non- LTE radiative transfer code RADEX 
(Ivan der Tak et al.ll2007h that computes the line intensities of a 
molecule for a set of physical parameters: kinetic temperature, 
H 2 gas density, molecular column density, background temper¬ 
ature and line width. We used the maps of CO 1-0, 2-1, *^CO 
1-0, C'^O 1-0, and CO 16-15 to perform minimization fit 
for a range of kinetic temperatures (30-700 K), column densi¬ 
ties (10*^ - 4x 10 '® cm“^) and volume H 2 densities (10^ - 5x 
10'** cm“^). In the area covered by the CO 16-15 map we have 
enough observational constraints to determine those 3 free fit¬ 
ting parameters. Eor all calculations we used a fixed line width 
value of 3.5 km s '. We adopted a cosmic background radiation 
field of 2 .73 K and used the molecular data from the LAMDA 
database (ISchoier et al.ll200^ . 

Eigure |9] shows the resulting distributions of the gas tem¬ 
perature, density, and column density of CO. As the CO 16-15 
line has an upper level energy of 750 K and a critical density 
of 10® cm~^ (see Table [TJ, one would expect that the strong line 
drives the solution to high temperatures and densities in the order 
of the critical density. However, it turns out that the low-7 lines 
can only be fitted with lower temperatures. To excite the CO 16- 
15 line, RADEX then has to increase the gas density to values 
in the order of 10* cm^*, a value much higher than derived in 
all other ways. This is a strong indication that the assumption 
of uniform gas parameters within the individual pixels of our 
map is not justified in S 140. It is much more likely that the true 
gas composition consists of a cooler component, mainly traced 
by the low-7 CO (and isotopologues) lines, while a second hot 
component is responsible for the CO 16-15 line. 

As all other lines discussed here are insensitive to gas densi¬ 
ties above 10® cm“*, the overestimate of the gas density obtained 
in this way should be irrelevant for the further analysis. The 
RADEX fit shows that the gas temperature peaks between IRS 1 
and IRS 2, roughly matching the CO 16-15 intensity map. The 
density structure shows a density peak at the positions of IRS 1 
and IRS 2. The column density follows mainly the distribution 
of the isotopologues which is higher towards the west region of 
the sources. It matches roughly the column density structure that 
was determi ned from EORCTAST, P ACS, and SCUBA dust ob¬ 
servations in iKoumnia et al.l (l2015l) . They also found a column 
density peak about 20" west of IRS 1 and IRS 2 and a smaller 
peak somewhat east of IRS 2. 

The CO column densities can be translated into abundances 
if we know the gas column from the dust observations of 
the central region. To derive gas column densities from the 
lOOt/m dust opacity map, we use the dust properties of model 
5 from lOss enkopf & Henning (1 19941) following the discussion 
in IKoumnia et al.l ( 20151) a ndthe standard con version factor of 
NuIAv = 1.9 X 10^' cm“^ (iBohlin et al.lll978l) . This provides a 
100 pm optical depth of 1.30 x 10“^'' cm^ x (Ah + 2 AH 2 ). The 
resulting column density map is shown in Eig. [10] 

The contours in the figure give the CO abundance com¬ 
puted by dividing the CO column from Eig|9]by the dust-based 
gas column density. Eor most of the cloud we find values of 
X(CO) = Aco/(A^h + 2 AH 2 ) around 2 - 3 x lO^”*. Lower values 
of about 1 X 10 are found at the submm-peak west of IRS 1. 
Somewhat lower values are also found towards IRS 2 and some¬ 
what higher values occur at the positions of IRS 1 and IRS 3. A 
reduction of the CO abundance towards IRS 2 would be consis¬ 
tent with the conversion of molecular material to atomic and ion¬ 
ized material seen in the bright fine-structure lines. The stronger 
decrease towards the submm peak could indicate CO freeze-out 
in very cold and dense material. But overall, the abundance vari- 
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Figure 10. Column density map of gas derived from FORCAST, 
PACS, and SCUBA continuum observations in the center of the 
region (colors), overlaid by contours of the CO abundance, run¬ 
ning from 1 X 10“"^ to 3 X 10“'*. 


ations within the map are about as large as the uncertainty from 
the dust columrQ. 

A CO abundance of more than 2 x lO '* would actually ex¬ 
ceed theJotal_amoum^ c arbon available in the gas phase (see 
e.g. ICardelli et al.lll99^ . iBlake et ^ (Il987h measured values 
X{CO) of only 6 x 10“^ towards the dense molecular material 
in the Orion molecular cloud and 2.5 - 3 x 10“^ in several clouds 
on larger scales, probably involving more atomic and ionized 
carbon. A higher CO abund ance of 1.3 x 10“ '* has been mea¬ 
sured towards NGC 2024 bv iLacv et alJ (Il994h . Our systemati¬ 
cally higher abundances can be explained by the uncertainty of 
the dust properties covering a fa ctor 2-3 in the 100 fim opac¬ 
ity (lOssenkopf & Henning|ll99^ . When assuming a dust model 
with a two times lower 100 jum opacity, the columns in Fig.fTOl 
increase by a factor two, the abundance values decr ease by the 
same f actor, and we obtain values close to those from iLacv et alJ 
(Il994ll . Unfortunately, the resolution of the low-7 CO maps is 
insufficient to determine the gas column density with the same 
resolution as the dust opacity. Hence, we can only conclude 
that the CO column densities are close to the upper edge of the 
range expected from the dust column densities. This means that 
most of the gas must be molecular so that we can use CO as a 
proxy for the total gas column in S140. This is consistent with 
the strong spatial confinement of the fine-structure line emis¬ 
sion in S140 indicating that CO provides a reasonable measure 
for the total gas column density, even if not all gas is molec¬ 
ular. Assuming X{CO) = 1.3 x 10“'* we obtain hydrogen col¬ 
umn densities Ah = Ah + 2Ah 2 between 7 x 10^^ cm“^ and 
1.97 X 10^^ cm“^ within the map covered by the CO 16-15 ob¬ 
servations. 

In the next step we can use the kinetic temperature and vol¬ 
ume density from the RAD EX fit of the CO lines (Fig.|9]) to com¬ 
pute the C^ and atomic oxygen column densities from the [C ii] 
and [Oi] observations. This allows to compare the abundance 
of the different species in the gas within the telescope beam 
along the line of sight. Actually, we do not expect that all species 
are spatially co-existing - we rather expect abundance gradients 
along the path and within the area of the beam - but the volume- 


^ The steep decrease of the dust column density and the correspond¬ 
ing virtual increase of the CO abundance at the boundaries of the map 
are probably d ue to missing flux at t he edges of the PACS footprint in 
the analysis bv lKoumnia et al](l201.5h . 
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Figure 11. Column density map of C^ derived from the mea¬ 
sured [Cii] intensities using the gas parameters from the 
RADEX fit to the CO lines. The contours show the [Cii] line 
optical depth a levels of t = 0.5,1.0,1.5,2.0. In the two pixels 
north of IRS 2, the optical depth is too large for a reliable column 
density estimate. 
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Figure 12. Column density map of atomic oxygen derived from 
the measured [O i] intensities using the gas parameters from the 
RADEX fit to the CO lines. The contours show the [Oi] line 
optical depth at levels of t = 0.05,0.10,0.15,0.20. 


integrated abundances can be easily compared with the results 
from more complex models and other Galactic and extragalac- 
tic observations. The CO 16-15 line provides strong constraints 
on the parameters derived from the RADEX fit and it is mainly 
stemming from hot and dense gas, the same gas that should also 
provide the main contribution to the fine-structure line emission. 
Therefore, the error made in assuming the same excitation con¬ 
ditions as seen by CO for deriving the C^ and oxygen column 
densities should be smaller than the factor two uncertainty that 
we have to face for the gas column anyway. To be consistent with 
the RADEX fit, we restrict the integration to the same velocity 
interval from -8.4 to -4.4 km s“*. Eig. [TTIand Eig. [TSlshow the 
resulting column densities and the optical depths of the lines. 
Although the line profiles of [O i] indicate optically thick line 
(Eig.lU, our resulting optical depths are r < 1. This is probably 
due to the strong self-absorption in the line creating an intensity 
dip for the considered velocity range, leading to too low column 
densities and optical depths. The [C ii] analysis results in col¬ 
umn densities 2.2x10*^ - 2x10*^ cm“^ and a range of optical 
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depths 0.28-2.31 with the higher values to be closer to IRS 2 in 
agreement with the estimate from the line profile in Sect. 13.11 

With the hydrogen column densities from CO the column 
density map translates into C^/H abundances, 2f(C^), between 
1.1 X 10“^ around the fine-structure line peak and 2 - 3 x 10“® 
for the rest of the map. For the oxygen abundance, we find values 
between 4 x 10“^ around IRS 2 and 1.0- 1.5 x 10“^ for the rest 
of the map. For both species the abundance at IRS 1 is as low as 
in the rest of the map outside of the fine-structure line emission 

peak. _ 

_ These values can be compared to PDR models (iRollig et alJ 

l2007h . At a PDR surf ace, almost all car bon is in the form of 
C^. For diffuse clouds. [Sofia et al.l (l2004h measured a abun¬ 
dance of 2f(C^) = 1.6 X 10“^, but in dense regions a larger frac¬ 
tion of carbon is incorporated in dust and PAHs leading to a 
PDR surface abundance o f 2f(C^) = 1.2 x 10 in dense clouds 
jWakelam & Herbs3l2008h . Deeper in the clouds where the UV 
is sufficiently shielded, turns into atomic carbon and sub¬ 
sequently into CO, so that the C^ abundance falls by three to 
four orders of magnitude, depending on the cosmic ray ioniza¬ 
tion rate. Even our value for IRS 2 falls much below the PDR 
surface limit, confirming that only in a small gas fraction car¬ 
bon is ionized through UV radiation while most of the gas is 
molecular. The column of molecular material should be at least 
ten times deeper than the C^ column produced in a PDR around 
IRS 2. 

Depending on the local conditions the oxygen abundance in 
dense clouds should not vary by more than a factor three. In the 
high-UV field limit, almost all gas-phase oxygen is in atomic 
oxygen, resulting in an abundance of X{0) = 3 x 10 relative 
to hydrogen. Deep in molecular clouds, most oxygen is incorpo¬ 
rated in CO so that only V(0) = 1 x 10 is left in the gas phase. 
Our measured abundance still falls below this lower limit. This 
can be explained by the fact that most oxygen is probably at 
temperatures below the range of 50-70 K fitted by RADEX for 
a homogeneous gas. The colder oxygen does not contribute to 
the fitted emission, but rather creates the absorption seen in our 
[Oi] line profiles. 

4.2. The [Cii] and [Oi] peak 

To estimate the size of the region responsible for the fine- 
structure line emission peak, we fitted the observed [O i] inte¬ 
grated intensity map by a Gaussian intensity distribution peak¬ 
ing close to IRS 2 convolved with the telescope beam at 4.7 THz. 
The [O i] map provides the main constraint to the source size due 
to the smaller beam compared to [C ii] and the CO lines. We al¬ 
lowed for a 2 arcsec pointing jitter in the fit when adjusting the 
location of the peak. The best fit leading to a relatively smooth 
background after the source subtraction is obtained for a source 
EWHM of 8.3" and a peak intensity of 76 K kms '. Eigure[T3] 
(top panel) shows the result. We see that the subtracted source 
covered 80% of the intensity at the position of the peak, but 
that there remains a somewhat more extended emission structure 
north of IRS 1. IRS 1 is part of that structure, but shows no con¬ 
centrated [Oi] emission associated with the infrared source. In 
contrast, we find a secondary [Oi] peak slightly north of IRS 3, 
the weakest infrared source in the field. 

A priori, the spatial extent of the emission in [O i] and [C ii] 
can be different, but if we subtract a source with the same size 
and the peak intensity of 212 K kms * from the [Cii] map 
we also obtain a very smooth map with no further indications 
of the IRS 2 contribution (Figure [13] central panel). This sug¬ 
gests that we see about the same gas responsible for the peak in 
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Figure 13. Intensity map of [O i] (top), [C ii] (center), and CO 
16-15 (bottom) after subtracting a Gaussian intensity distribu¬ 
tion of FWHM= 8.3" from the observed map. The overlaid con¬ 
tours show the original intensity maps at levels of 8, 24, 40, 56, 
72 K kms-* ([Oi]), 50,70,90,120,150,180K kms-' (CII), and 
10, 20, 30,40, 50, 60 K (CO 16-15). 


[O i] and [C ii] spatially well confined to a region close to IRS 2. 
Repeating the same experiment for the CO 16-15 map, provides 
a peak intensity of 46 K km s ', but a clear residual of emission 
from IRS 1 and IRS 3 (Figu re [T3I bottom panel). This agrees 
with the previous analysis of iKoumnia et al.l (12015h finding ex¬ 
tended hot molecular gas around IRS 1. 

The identification of a spatially well confined emission struc¬ 
ture close to IRS 2, responsible for the observed strong fine- 
structure line emission allows to quantify the properties of the 
emitting gas. With an FWHM of 8.3" we obtain a spatial angle 
of the emitting region of 1.8 x 10-® sr and a diameter of 0.03 pc. 
This allows to compute absolute luminosities from the integrated 
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intensities and the source projected area of 1.0 x 10^^ cm^ as¬ 
suming isotropic radiation. For the [Oi] line we obtain 0.28 L© 
and for the [C ii] line 0.05 L© . This can be compare d to the total 
energy input at that location. iKoumnia et alJ (1201 5h estimated a 
total luminosity of the source at IRS 2 of 2000 L©. This means 
that only 0.016% of the source luminosity is radiated away in 
the two fine-structure cooling lines. We will further discuss this 
in Sect. 14.41 

For the CO 16-15 line we obtain a corresponding luminos¬ 
ity of only 0.01 L©, but over the full CO ladder this will add up 
to a considerable net cooling. If we make the extremely conser¬ 
vative assumption that all CO line intensities are constant at the 
value of the 16-15 transition for lower J transitions and zero for 
higher transitions, we obtain a lower limit for the CO cooling 
luminosity of 0.045 L© already matching the [C ii] luminosity. If 
we use the full SED ladder fit to the narrowest line component 
from the next section (14.31) . we obtain instead a cooling power of 
0.096 L© and if we include the emission over the total line width, 
the CO luminosity sums up to 0.13 L©, i.e. a value that is about 
half the cooling power in the [Oi] line. In spite of the prominent 
fine-structure line peak at IRS 2, we thus find a situation where 
the molecular line cooling is relatively important next to the fine 
structure line co oling. This is simil ar to the situation found in 
other protostars (iKarska et al.ll2013l) . Obviously, the situation is 
even more extreme for IRSl, where we measure strong high-7 
CO emission, but basically no [C ii]. 

We can estimate the gas column density from the integrated 
[C ii] line intensity, by using the relativ ely constant rCn] emis - 
sivity per column density (Eq. (2) from lOssenkopf et akl 2013h . 
Eor a temperature of about 200 K we obtain: 




edv ^ 560 


Kkms ' 
cm^^ pc 


xA^c+ 


( 1 ) 


As the emissivity is only weakly temperature-dependent, this 
value applies within a 15 % error bar to all excitation temper¬ 
atures between 100 K and 500 K, covering the typical range ex¬ 
pected in PDRs. Eq. ([1]) provides a direct translation of intensity 
to column density for optically thin emission. Because we know 
the optical depth from Sect. 13.11 the equation can also be used 
with the corresponding optical depth correction of 45 % to com¬ 
pute the C^ column density for IRS 2. 

We obtain a C^ column density of 1.7 x lO'^ cm“^. This 
is slightly lower than the value of 2 x 10'^ cm“^ obtained for 
the total [C ii] emission towards IRS 2 measured in the RADEX 
simulation for the full map in Sect. 14.11 However, here we con¬ 
sider only the 73% of the [Cii] emission at that position that 
are attributed to the compact source ignoring a small part of 
the total column density. Using t he carbon gas-phase abund ance 
of 1.2 X 10 for dense clouds dWakelam & Herb^l2008ll (see 
Sect l4.1l) and assuming that all carbon is ionized, the column 
density of 1.7x 10'^ cm“^ provides a lower limit for the total col¬ 
umn density of the [C ii] emitting gas of 1.4 x 10^^ cm“^. Higher 
values result for partial ionization. If we assume that the [Cii]- 
emitting gas fills the observed emitting volume of 1.0x10^' cm^, 
we deduce a minimum gas density of 1.4 x 10^ cm“^. The de¬ 
rived optical depth of the clu mp seen in dust emis sion close to 
IRS 2 is, however. Ay ~ 100 dKoumnia et al.ll2()T^ correspond¬ 
ing to gas column densities of 2 x 10^^ cm“^. That means, that 
only 7 % of the gas is emitting in [C ii]. This could be explained 
by dumpiness of the gas within the 8.3" peak and high-density 
clumps, or more naturally, by a PDR layering structure where 
the [C ii] emission only stems from gas with Ay S 2 facing the 
embedded sources. 


Table 3. Physical conditions as derived with RADEX for se¬ 
lected positions 


Position 

Tkin 

[K] 

N(CO) 
[10“* cm-2] 

Nh 

[10^- cm-2] 

«H2 

[cm“^] 

IRS 1 

75 

6.6 

11 

3x10*' 

IRS 2 

73 

16 

27 

7.1x10* 

Interface 

40 

1.0 

1.5 

4x10^ 


Notes. The synthetic integrated intensities are calculated using a 
FWHM of 3.2 km s^* for IRS 1 & 2.7 km s^* for IRS 2. The parameters 
towards the IF were derived in addition the CO 1-0 & CO 2-1 . 
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Figure 14. CO ladder from the RADEX fit to the line intensities 
at the three prominent positions. 



Unfortunately, it is not possible to perform the same estimate 
for the [Oi] emission. The high-temperature LTE limit of the 
[Oi] emissivity is about three times as large as for [Cii] (Eq.[TJ, 
but it only applies to densities above 10® cm^^ and temperatures 
above 200 K, and it also requires optically thin emission or an 
accurate estimate. The line profiles showed already that [Oi] is 
optically very thick and the intensity of the [Oi] line is three 
times lower than that of [Cii]. To check the consistency of the 
parameters in the frame of a real radiative transfer computation 
we perform RADEX runs in the next section. 

4.3. RADEX fit of individual sources 

A more accurate estimate can be obtained when using the known 
gas parameters from a fit to the observed CO transitions and tak¬ 
ing the line optical depth into account in the frame of a RADEX 
computation. As all lines are bright towards the central sources, 
IRS 1 (0,0) and IRS 2 (0,20), we can avoid the low-J CO lines 
there that show complex line profiles and are very optically thick 
and fit only the higher transitions of CO and all available transi¬ 
tions of the isotopologues in RADEX. Eor the int erface (-40,54), 
the RA DEX analysis was already performed in iKoumnia et al.l 
(1201 5h using the low-7 CO lines and an additional HIEI cut mea¬ 
suring CO 9-8 and '^CO 10-9. As our CO 16-15 observations do 
not cover this region, no further constraints are provided to the fit 
so that we can reuse the results obtained there. The HIEI cut also 
covered IRS 1 so that we have the most complete data set for 
this source. Table |2] contains all transitions and line properties 
that were used for the RADEX fit and the further analysis. 

Eigure [14] shows the resulting CO SED from the RADEX 
fit for all lines up to 7 = 17-16 and Table [3] summarizes 
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Table 2. Line parameters of the narrow component of CO and isotopologues towards IRS 1 and IRS 2 - GREAT, IRAM & HIFI 
data. 



IRS 1 

IRS 2 

interface 

Line 

Vlsr 

FWHM 

J T„i,4v 

Vlsr 

FWHM 

J T,„i,4v 

Vlsr 

FWHM 

J T„,i,dv 


[km s“’] 

[km s“*] 

[K kms“‘] 

[km s“’] 

[km s“*] 

[K kms^*] 

[km s“'] 

[km s“*] 

[K kms^*] 

CO 9-8 

-5.78 ±0.01 

4.87 ± 0.06 

117±2 




-7.80 ± 0.07 

2.50 ± 0.20 

1.20 ±0.20 

CO 13-12 




-6.76 ± 0.03 

3.30 ± 0.08 

71.1 ± 1.4 




CO 16-15 

-6.55 ± 0.06 

3.5 ±0.2 

40.3 ± 1.7 

-6.43 ± 0.04 

2.7 ±0.1 

37.2 ± 1.3 




‘^CO 1-0 

-6.45 ± 0.03 

2.48 ± 0.09 

28.7 ± 1.9 

-6.54 ± 0.06 

2.6 ±0.1 

49.4 ± 1.8 

-8.12 ±0.37 

1.34 ±0.70 

0.84 ± 0.26 

‘*CO 10-9 

-6.71 ±0.02 

2.41 ± 0.08 

12.8 ± 1.1 






< 3cr 

C'*0 1-0 

-6.52 ± 0.05 

2.4 ±0.1 

6.8 ±0.6 

-6.64 ± 0.04 

2.36 ±0.1 

8.4 ±0.1 



< 3o- 

C"*0 9-8 

-6.7 ±0.1 

2.6 ±0.2 

2.8 ±0.3 






< 3cr 

[Cii] *P3/2-^ Pm 

-7.44 ± 0.04 

5.94 ± 0.09 

72 ± 1 

-6.27 ±0.01 

4.62 ± 0.03 

159.5 ±0.9 

-7.92 ±0.01 

3.24 ± 0.02 

119.6 ±0.7 

[Ol] V, -3P2 

-5.9 ±0.2 

5.3 ±0.4 

9.5 ±0.7 

-6.4 ±0.1 

5.9 ±0.2 

25.4 ± 0.7 

-9.1 ±0.6 

2.7 ±0.5 

3.1 ±0.9 


Notes. The line parameters are measured in a 21" Gaussian aperture around the indicated position (telescope beam + convolution kernel) to match 
the coarse resolution of the *^CO 1-0 map. 


the resulting gas parameters. In the fourth column, we trans¬ 
lated the CO column density into an equivalent hydrogen column 
density using the abundance factor of 6 x 10“^ from Sect l4.ll 
The higher column density towards IRS 2 compared to that of 
IRS 1 is driven by the spatial distribution of the low-J CO iso¬ 
topologues and the high-J CO that is stronger towards IRS 2 
(see Fig. |6li. The hydrogen column densities towards the two 
sources approximately match the visual extinction of Ay ~ 50 
towards IRS 1 , and Ay a; 100 towards IRS 2 measured by 
iKoumnia et al.l (l2015h when assu ming a standard con version fac¬ 
tor of NnjAv = 1.9 X 10^' cm“^ (iBohlin et al.lll978h . 

The IRS 1 and IRS 2 are warmer than the interface and have 
a similar temperature. The absolute value of the temperature is, 
however, questionable as it depends on the assumption of the 
uniform gas parameters within the beam. As discussed already 
in Sect. 14. II this also leads to very high densities. The volume 
density estima tes for IRS 1 and IRS 2 are much higher than the 
ones found in iKoumnia et al.l (l20L5h (~ 10 ® cm“^) due to the 
inclusion the high-J CO transitions in the fit (13-12, 16-15). 
We should rather expect that the different transitions arise from 
different layers around the infrared sources. The fit produces a 
good match to all observed lines, so that we have an accurate 
estimate of the total cooling through all CO lines from the ob¬ 
served sources. 

We use again the RADEX gas parameters to compute the 
C^ and oxygen column density from the line intensities given 
in Table |2 For IRS 1, we get column densities for C^ and O 
of 7 X 10'^ cm“^ and 1.8 x lO'^ cm“^, respectively, for IRS 2 
5 X 10** cm“^ cm“^ and 6 x 10** cm“^ cm“^. When comparing 
these values with the map analysis, we notice that the C^ col¬ 
umn density towards IRS 2 is increased by a factor two because 
we use the full line width now and extract the intensity exactly 
at the position of the fine-structure line peak. All other values 
are basically unchanged. At the interface, we obtain a C^ col¬ 
umn density of 5 x 10** cm“^ and an oxygen column density of 
3 X 10*® cm“^. We can compare these values with the column 
density that we obtain from the total gas column density mea¬ 
sured by CO and the upper limits of the abundance of the two 
species discussed in Sect l4.ll In this way we obtain column den¬ 
sities of A(C^) = 2 X 10** cm“^ and A^(0) = 5 x 10** cm“^, i.e. 
much lower than computed here. The explanation must be given 
by the PDR nature of the cloud interface. As most CO is disso¬ 
ciated there, it is no longer a good measure for the total column 


density. The actual column density must be at least a factor six 
higher when using the gas temperature obtained from CO. 

4.4. Cooling balance 

The ratio of the energy emitted by the interstellar gas in the two 
main far-infrared cooling lines ([O i] and [C n]) relative to the in¬ 
tegrated far-infrared continuum emitted by the dust can be used 
to estimate the gas-heating efficiency when assuming that the 
far-infrared flux reflects the total energy of the incident stellar ra¬ 
diation held by converting it to infrared-wavelengths through ab- 
sorption and re-em ission in an optically thick medium (see e.g. 
lOkada et al.ll2013|) . As access to the [Oi] line became possible 
only very recently, the [C ii]/FIR intensity ratio is often used to 
measure the efficiency of the gas heating. This allows to address 
variations in the efficiency of the photoelectric heating from dust 
grains because photoelec tric heating is the most imp ortant gas 
heating process in PDRs (iHollenbach & Tielenslll99^ . 

Although the definition of the far-infrared varies be¬ 
tween different papers (e.g. as FIR integrated from 42 
to 122 jum, or as TIR integrated from 3 to 1100 gm; 
[Dale & Heloul l2002l) . the typical range of the gas heat¬ 
ing efficiency in Galactic PDRs, star-forming regions in 
EMC and M33 is 10~*-10~^ both i n case of being 


traced by [Oil+rCiil (Okadaetal. 20131 

Lebouteiller et al.l 

l2012t iMookeriea et alJ 201 ll Mizutani et al. 

2004h and 

traced 

only by ICiil (lOkada et al.l 12015; Lebouteiller et al. 

2012); 

IMookeriea et alj|201 ll). As an extreme case.lVastel et al. 

(200ll) 

found an efficiency of 10 "* in W49N, which is illuminated by an 


intense UV held. 

For the region around the embedded sources, the far-infrared 
luminosity w as derived from PACS a nd SOFIA-FORECAST ob¬ 
servations bv iKoumpia et al.l (l2015h . It was integrated from the 
observed intensities between 11 and 187 gm further extrapo¬ 
lated linearly to zero flux levels at longer wavelengths. FigurefTSl 
compares the energy radiated in our fine structure lines with the 
total infrared flux (d > 11 fim) at the resolution of 13" from 
the PACS continuum dat^^ In the background we see the in¬ 
frared flux from the three embedded sources with luminosities 
of 10000 L ^dRS 1), 200()L^(I RS 2), and 1300 Lq (IRS 3), re¬ 
spectively (IKoumpia et al.ll20f^ . The contours show the ratio of 

^ Our line cooling estimates for the 8.3" source at IRS 2 in Sect. 14.21 
slightly deviate from the numbers given here because of the coarser 
resolution applied here compared to the direct fit of the emission profile. 
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Ac( [orcsec] 

Figure 15. Comparison of the line cooling with the infrared 
continuum cooling. The background colors show the total in- 
frared flux integrate d from wavelengths starting at 11 fim 
(iKoumnia et al.ll2015l) . The contours give the decadic logarithm 
of the ratio of the integrated line flux relative to the continuum 
flux. In the top panel we use the sum of the [O i] and [C ii] fluxes, 
in the lower panel, we only show the [C n]/TIR ratio as used in 
many existing investigations. 

the sum of the [O i] and [C ii] fluxes relative to the total infrared 
(upper plot) and the [C ii] line only (lower plot) in logarithmic 
units. 

The [Cii]/TIR ratio shows an almost spherically symmet¬ 
ric picture around IRS 1, being dominated by the variation of 
the infrared continuum flux. The values of 10“^-10“^, typically 
observed in other Galactic PDRs, are only seen at the south¬ 
western edge of the mapped area closest to the illumination 
by HD 211880. Around the embedded sources, the infrared lu¬ 
minosity grows without a corresponding growth of the fine- 
structure line luminosity resulting in [C ii]/TIR ratios below 10“^ 
close to IRS 1. The fine structure peak at IRS 2 only weakens this 
global scenario slightly, resulting in a “dip” in the contours next 
to IRS 2. 

When including the [Oi] line in the comparison, the ratio 
increases by a factor three close to the sources and the fine- 
structure line peak close to IRS 2 creates a plateau in the en¬ 
ergy ratio a about 2 x 10“"^ in agreement to the results ob¬ 
tained in Sect l4.2l At lower densities towards the map bound¬ 
aries the [Oi] contribution is negligible. For IRS 1, even the 
cooling through the lines of the CO ladder is four times stronger 
than the sum of the [C ii] and [O i] cooling power (see. Fig. [T4l) . 
for IRS 2 CO sums up to half of the fine-structure line cool¬ 
ing. Including the CO rotational ladder in the gas cooling bud¬ 
get raises the line-to-continuum ratio to 1 - 3 x lO^'*, i.e. val- 



60 40 20 0 -20 -40 -60 -80 

Aa [arcsec] 


Figure 16. Comparison of the total line cooling flux provided 
by the sum of the [O i] and [C ii] lines (colors) to the [O i]/[C ii] 
line flux ratio (contours from 0.4 to 2.8 in steps of 0.4). The 
computation is performed in a common effective beam of 15". 


ues that are still extremely small in Galactic context and that 
would be interprete d as FIR line deficit when observed in ex¬ 
ternal galaxies (e.g. Luhrnan et al.iri998l: iMalhotra et al.l l2()0Tt 
iHerrera-Camus et^ l2015h . However, none of the proposed ex¬ 
planations for the FIR line deficit in external galaxies applies to 
our case (see Sec.|5]l. 

The ratio between the [O i] and [C ii] intensities is considered 
as one of the main characteristics of PDRs. For high UV fields, 
the gas temperature is always above the upper level energy of 
both transitions so th at the ratio between t he two lines directly 
measures the density (iKaufman et al.l 199^ . However, due to the 
higher energy of the [O i] transition it falls off quite fast at low 
UV fields resu lting in negligible [O i] cooling for low densities 
and UV fields (iRollig et al.ll2006h . 

As visible already from Fig. the relative contribution 
of [O i] and [C ii] to the cooling changes significantly over the 
cloud. Figure [16] shows the total line cooling budget of the [Oi] 
and [C It] lines in colors and the [O i]/[C ii] line ratio in contours. 
We And a relatively constant ratio of 2.5-3 in the whole area 
of the central cluster connecting IRS 1 and IRS 2. Here, [Oi] is 
dominating the gas cooling. A value of three is representative for 
either very high UV fields and low densities or low UV fields and 
very high densities. For the high densities and UV fields around 
the central cluster, we would however expect much higher val¬ 
ues due to the high abundance and emissivity of oxygen in that 
regime. 

When leaving the area of the central cluster, the contribu¬ 
tion of [O i] to the total line cooling budget quickly drops off. 
This should reflect a combination of a reduced density and tem¬ 
perature as measured by CO (Fig. |9]l. While the density drops 
more quickly south of IRS 1, the temperature drops more quickly 
north of IRS 2. The asymmetry of the [O i]/[C ii] ratio plot indi¬ 
cates that the density is the dominating factor here. In the bulk 
of the cloud and the interface, the cooling strength of the [O i] 
line can be completely neglected. [C ii] is the only relevant gas 
cooling line there. 

4.5. PDR model interpretation 

To overcome the RADEX assumption of constant excitation 
conditions within the emitting region as an oversimplification 
we switch to theoretical models for PDR structures that self- 
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consistently compute the temperature and chemical structure of 
PDR layers around illuminating sources, allowing us to com¬ 
pare the line observations with the model predictions for differ¬ 
ent densities and radiation fields. 

One critical input parameter is the strength of the imping¬ 
ing UV radiation field. It can be computed by assuming that in 
an optically thick configuration all UV photons are eventually 
converted into observable infrared photons. Then we obtain the 
strength of the incident radiation field in terms of the Draine field 
(lDrainelll978h as: 


^ 2 X 2.7 X lO-^W/m^ 

where the factor two in the denominator accounts for half of 
the heating photo ns at wavelengths out of the UV band (e.g. 
[Rollig etdl I 2 OI ill . Using the total infrared flux values from 
Sect. 14.41 results inx = 5.6 x IQ'^xo for 2 and 1.2 x 10^;^^o 
for IRS 1. For the interface, the UV field has been derived pre¬ 
viously to_^ll_between_J_402P([Poei^n[&^Spaa^|2005l) and 
llTimmermann et aOl996h. 

As a simple and straight-forward model, we can treat 
the PDRs in S 140 as plane-parallel slabs simulated by 
iKaufman et al.l (Il999t) . In the model line and continuum inten¬ 
sities are self-consistently computed for a face-on configuration 
taking the optical depths of the lines into account. This allows 
us to directly read the source parameters from the [C ii] line in¬ 
tensity, the [Oil/[C III intensity ratio, and the ([Oi]-i-[Cii])/TIR 
ratio (Figs. 3-6 from IKaufman et ^ 1 19991) . As the model has 
only the two parameters gas density and incident radiation field, 
the problem is over-determined allowing us to assess the consis¬ 
tency of the description. Depending on the source geometry, we 
may, however, introduce the inclination angle as an additional 
free parameter providing a line-of-sight increase of the column 
densities because the model is only computed for a face-on PDR. 

This model is clearly appropriate for the external cloud 
surface representing an almost edge-on PDR illuminated by 
HD 211880. The main [Cii] and [Oi] emission from a plane- 
parallel PDR stems from gas within Ay < 2 from the cloud 
surface (iRollig et alJl2007l) . This is also the approximate scale 
of ch emical stratification in the PDRs (iHollenbach & TielensI 
Il999l) . In Fig. |2] we can directly measure the stratification be¬ 
tween CO 2-1 and [Cii] as 15" or 0.05 pc. This allows to 
estimate the gas density as hh ~ 2 x (Vh/Av)/0.05 pc = 
2.5 10^ cm-\ 

When reading the model predictions for the given parameter 
combination of density and previously fitted UV field for the in¬ 
terface, we obtain a [Cii] line intensity of 3 x 10'^ W sr“* nT^, 
a [Oi]/[Cii] energy ratio of 5, and a line-to-continuum ratio 
of 0.02. Unfortunately, we do not have numbers for the in¬ 
frared continuum flux at the interface position as it was out¬ 
side of the region observed by FORECAST and PACS. We can, 
however, assume that the ([Oi]h-[Cii])/TIR ratio falls above the 
value of 10“^ seen in the south-western edge of the mapped 
area (Fig. dsii. The measured [Cii] line intensity is, however, 
10 ® W sr“^ m“^, a value that is obtained in the model only for 
UV fields above 2 x 10‘^xo and high densities. In contrast, the 
measured [Oi](63 jum)/[Cii] energy ratio is 0.3, a value that is 
found in the model only for UV fields of about 30;^^o and below. 

A solution can be achieved by introducing the inclination 
angle of the PDR as additional parameter and changing the as¬ 
sumption on the impinging UV held. If the PDR inclination pro¬ 
vides a geometrical magnification of the [C ii] intensity in the 
optically thin case and we assume that [Oi] is optically thick 
so that it is not amplified in the same way, the observed [C ii] 


intensity can be corrected downwards and the [Oi]/[Cii] ratio 
upwards. For an amplification factor of eight, we obtain a match 
of both quantities with the model at a UV held of 60 xo- This 
is lower than obtained in the previous PDR model fits, but still 
six time s highe r than measured deeper in the S 140 region by 
iLi et al.l (l2002h . Because our gas density value is based on a 
rough estimate only, we should assume a a factor two error bar 
for all numbers for a conservative approach. The geometrical 
amplification of a factor eight is obtained by an inclination angle 
of 80-85 degrees relative to the face-on orientation. The study of 
the [Oi](145 jum)/[Oi](63 yum) line ratio based on PACS data in 
Appx.lBlsuggests a somewhat lower amplification factor of four. 

We can try to use the same model as an approximation for the 
internal PDR around the cluster at IRS 2 that gives rise to the ob¬ 
served [C ii] and [O i] peak. Here, we assume that the PDR will 
form a shell around the cluster with a face-on front and back side. 
For the gas density, the value of 1.4 x 10^ cm“^ from Sect. 14.21 
provides a lower limit. Reading the observables from the model 
for this density and the UV held of 5.6 x IQ^xo gives a [C ii] line 
intensity of 1 x 10“® W sr“' m“^, a [O i]/[C ii] energy ratio of 40, 
and a ([O i]-i-[C ii])/TIR ratio of 3 x 10^^. Comparing this to our 
measured values of 1.5 x 10“® W sr“'m"^, 5.6, and 2 x 10“"^ we 
And a good agreement for the [C ii] line intensity when assum¬ 
ing that the radiation from both surfaces adds up in the optically 
moderately thick line (f = 2, see Sect. 13.11) . However, the ob¬ 
served [O i] line intensity and is about seven times lower than 
predicted resulting in a line-to-continuum ratio that is also six 
times lower. 

The much lower [O i] intensity in the observations may stem 
from the shell configuration where we expect a decreasing [O i] 
excitation temperature towards the observer in contrast to the 
face-on configuration of the PDR model where we directly ob¬ 
serve the illuminated side of the PDR. For [C ii] both configu¬ 
rations do not make a big difference as C^ is only abundant in 
the narrow Ay < 2 layer, but atomic oxygen is also abundant 
in the cool environment at the back side of the PDR, potentially 
absorbing a large fraction of the [O i] emission from the PDR. If 
we assume that this self-absorption is as strong as a factor seven, 
much higher than estimated from the line profiles in Sect. 13.11 
we obtain a self-consistent pictur e for the PDR of IRS 2 in 
the frame of the PDR model from IKaufman et aP (Il999h . This 
self-absorption by a factor of seven is also consistent with the 
[Oi](145 jum)/[Oi](63 jum) ratio observed by PACS (Appx.IBI). 

By further inspecting the model output as a function of den¬ 
sity, we can also constrain the maximum density of the PDR. 
When increasing the density by a factor ten above our lower 
limit, the [Cii] intensity drops by a factor two while the pre¬ 
dicted [Oi] intensity remains approximately constant. In the 
shell configuration this cannot be compensated by a variable in¬ 
clination angle so that we can constrain the actual density range 
to densities below 10® cm“^. This is consistent with the cen- 
tr al density estimate of 9 x 10® cm ® from the dust modelling 
in iKoumpia et al.l (l2015h . With the assumption of a very strong 
[O i] self-absorption, the fine-structure line peak at IRS 2 is thus 
still consistent with the standard PDR picture in a UV held of 
5.6 X lO^xo and a density between 1.4 x 10® and 1 x 10® cm“®. 

However, we do not succeed any more in providing a con¬ 
sistent PDR model explanation for the fine structure line deficit 
from IRS 1. The combination of the [Cii] intensity of 4 x 
10“® W sr“* the [Oi]/[Cii] energy ratio of 1.6, the line- 
to-continuum ratio of 2 X 10 ®, and the known UV held of 
1.2 X 10®;t'o provides a solution for a gas density of 50 cm“® 
and a factor 1.5 self-absorption of the [Oi] 63 /tm line only, but 
this low density is in clear contradiction to the analysis of the CO 
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lines i n Sect. 14. II and the dust emission profile in lKoumnia et alJ 
(l2015h . In fact, the density around IRS 1 should be at least as 
high as towards IRS 2. 

5. Discussion 

Our observational results raise two main questions: 

- Why is the fine-structure line cooling from IRS 1 so weak 
compared to the continuum, with [C ii]/TIR ratios as low as 
10 -^? 

- Why do we see significantly more fine-structure line emis¬ 
sion from the weaker sources IRS 2 compared to IRS 1 ? 

We can discard a number of easy explanations for the lack of 
[C ii] emission from IRS 1 and the relative weakness of the [O i] 
and [Cii] emission from IRS 2: 

High resolution radio continuum observation s using the 
VLA and MERLIN (iTofani et aI.lir^lHoarell200^ showed ul¬ 
tracompact H ii regions at the positions of IRS 1 and IRS 2. This 
shows that the embedded sources produce enough UV photons 
with energies above 13.6 eV to ionize hydrogen. Therefore, they 
should also be able to ionize carbon (11.2 eV) and PAHs and 
dust grains leading to production and a sig nificant photo¬ 
electri c heating of the gas. In the same direction iMalhotra et al.l 
(1200 Ih proposed radiation from an older stellar population with 
an overall redder spectrum to explain the [C ii] deficit, but in our 
sources we can be sure to find a rather very young cluster. 

PACS observations of the [N in] 57 /rm line in the f rame of 
the Herschel key project WADI dOssenkonf et al.ll201 ih did not 
provide any detection of the line in 47" x 47" area centered at 
IRS 1. Because of the similarity of the second ionization poten¬ 
tial of nitrogen and carbon, the lack of rules out that car¬ 
bon is also in the form of not contributing to the [C ii] line. 
Together with the non-detection of the [O iii] line, this rules out 
that a high ionization parameter explains the lack of [C ii] and 
[O i] emission. Unfortunately, the existing Spitzer IRS spectra do 
not allow to constrain the abundance of other ionized species as 
they are saturated in most of the field due to the large continuum 
flux. Future observations of infrared lines might show whether 
there is a general line deficit from the Hit regions that might 
be due to strong UV absorption by dust within the H ii regions 
preventing any FUV photons from arriving at the molecular gas 
around them. 

The intensity of the far-infrared lines may be reduced 
by dust extinction for very large column densities (see e.g. 
lEtxaluze et al.ll2()T^ . Mai nly constrained by 450 jjm SCUBA 
observations we derived in iKoumnia et al.l (120 15h a total visual 
extinction of Ay ~ 50 towards IRS 1, and of up to Ay 100 
towards IRS 2. The dust opacities should fall between the values 
for diffuse clouds (82 cm^/g at 63 //m, 12 cm^/g at 158 jum, 
iLi & Draind 1200 ll) and those fo r dense clumps (211 cm^/j ' 
at 63 jjm, 41 cm^/g at 158 um. lOssenkonf & Henning! 1 1994 . 
model 5) providing optical depths of T 63 ^nj = 0.1 - 0.25 and 
Ti 58 ^m = 0.01 - 0.05 for IRS 1 and T 63 ;jm = 0.2 - 0.5 and 
Ti 58 ^m = 0.03 - 0.2 for IRS 2 when we assume that half of the 
dust is located in front of the line-emitting gas. This may provide 
a reduction of up to 40 % for the [O i] intensity towards IRS 2, 
but only a negligible effect for all other cases, in particular for 
the rC III line, not explaining the overall large deficiency. 

[Herrera-Camus et^ (1201 5h proposed a redistribution of the 
cooling power from [C ii] to the [O i] 63 fim line in dense ma¬ 
terial as an explanation for the [C ii] deficit, but as we observed 
both lines, we can be sure that this can only apply for the self- 
absorbed part of the [O i] line as all other photons are measured. 


We find, however, a redistribution of cooling power to the CO 
rotational lines. They can carry as much as half of the energy of 
the [Oi] line, but clearly do not account for a factor 100 differ¬ 
ence. 

Based on the spatial and spectral information in our data 
we can exclude a shock-origin of our fine structure peak. It is 
about 5" apart from the outflow knots from IRS 1 mapped by 
I Wei gelt et al.l (l2002h and the fine-s tructure lines do not s how the 
outflow wings seen e.g. in [C i] bv iMinchin et al.l (Il994t) or seen 
in our low-7 CO data. In case of outflow-excited [O i] emission 
we expect a wider line width in [Oi] compared to CO, as ob¬ 
served e.g. in G5.89-0.39 by iLeurini et al.l (1201 5h and not the 
opposite behavior. Moreover, the [O i] line has a slightly blue- 
skewed profile that is rather representative for infall than for out¬ 
flow. 

In the frame of PDRs, low line-to-continuum cooling ratios 
are explained by either very high or yery low charging param¬ 
eters ;^f/n (e.g. iKaufman et allll999l) . High ;^f/n values lead to 
positively charged dust grains and PAHs, lowering the photo¬ 
electric efficiency . They are obtained for very high UV fields 
at low densities. In the opposite case all material remains cold 
and no C^ is produced so that most energy is emitted in CO 
lines. However, our CO and dust observations exclude low den¬ 
sities around the infrared sources (Sect. l43T l. the measured [C ii] 
intensity excludes the cold-gas solution, and the UV fields are 
constrained by the measured infrared luminosity. 

Alternatively, the photoelectric efficiency can be reduced by 
the destruction of PAHs s ince PAHs are the dom inant carrier of 
the photoelectric heating (iBakes & Tielensiri994ll . As mentioned 
above, the saturation of the IRS spectra around infrared sources 
prevent us to examine the PAH spectra directly, but when we 
compare the relative strength of the PAH 11.3 fim feature against 
the underlying continuum southwest of IRS 1, we find some re¬ 
duction in the area closer area to IRS 1 (20-30" from IRS 1) 
compared to 30-50" from IRS 1. However, the spectra show still 
an infrared excess from very small dust grains indicating that 
the reduction in the PAH abundance can only weakly reduce the 
overall surface for the photoelectric effect. 

The relative weakness of the fine structure lines and the low 
[O i]/[C ii] ratio must be explained by deviations of the condi¬ 
tions in S 140 from that “standard PDR”. For the [Oi] line the 
low intensity can be explained by a scenario with a decreasing 
[O i] excitation temperature towards the observer, resulting from 
a decreasing gas temperature and/or density. Due to the high 
column density of atomic oxygen at low temperatures, this ge¬ 
ometry can provide a strong self-absorption, explaining the fac¬ 
tor seven reduction of the [O i] line intensity observed towards 
IRS 2. However, this mechanism can hardly be responsible for 
the same kind of reduction in the [C ii] line intensity towards 
IRS 1. As quickly recombines in a dense gas at low tem¬ 
peratures, finally forming CO, we do not expect large quantities 
of cold in the line of sight that can produce the same kind 
of self-absorption for the [C ii] fine-structure line. A very spe¬ 
cial configuration would be needed that allows for a significant 
amount of cold low-density gas in a clumpy medium behind 
a layer of dense and warm gas facing the H ii region around the 
embedded cluster. 

Beyond the relative weakness of the fine-structure lines rel¬ 
ative to the continuum flux, we find puzzling combinations of 
detections: 

1. IRS 2, with 2000 Lq shows strong [Oi], [Cii], and CO (and 

isotopologues) emission at about -6.5 kms^' and has asso¬ 
ciated ultracompact H ii regions. 


14 














































Ossenkopf et al.: Fine structure line deficit in S 140 


2. IRS 1, with 10000 L© shows only weak [C ii] and [O i] emis¬ 
sion at the main velocity of the cloud of -8 km s *. The CO 
lines have a somewhat smaller intensity than on IRS 2 and 
the optically thin lines are centered at the IRS 2 velocity, 
showing only an emission wing at -8 kms indicating also a 
relative lack of CO emission from IRS 1. An ultracompact 
H II region and an outflow were detected. 

3. IRS 3, with 1300 L© shows [Oi], but hardly any [Cii], CO is 
comparable to IRS 1, but spatially not separated in our beam. 
The source, however shows no associated H ii region. 

6. Summary 

Our GREAT observations of S 140 showed a pronounced peak 
of the emission of the [O i] 63 yum and [C ii] lines 20" north of 
IRS 1, the main embedded source dominating the infrared field 
in the whole region. The fine-structure-line peak can be asso¬ 
ciated with a weaker infrared source, IRS 2, and the GREAT 
observations resolve the size of the emitting region to have a 
diameter of 0.03 pc. The velocity of the gas at that position is 
offset by about 1.5 kms^' from the bulk of the molecular ma¬ 
terial in S 140. The gas density of the emitting gas must be at 
least 1.4 X 10^ cm“^ to allow for the observed strength of the 
[C ii] line. In spite of its absolute strength, the [O i] 63 yum line, 
is however, relatively weak when compared with [C ii] and the 
infrared continuum. This can only be explained by a strong self¬ 
absorption in the foreground medium with a large optical depth. 
Eor [C ii] we measured an optical depth of t » 2. 

Erom IRS 1 we find only relatively weak [C ii] and almost 
no [Oi] emission in spite of its high luminosity of lO^L© and a 
know ultracompact H i region. The main heating source of dust 
and large scale-gas in S 140 is not responsible for the [Cii] and 
[O i] peak, but shows an extreme fine-structure line deficit. While 
for IRS 2, the low value line of the ([Cii]h-[Oi])/TIR ratio, also 
matching the criterion for a line deficit, can be explained by the 
strong [O i] self absorption, this explanation does not work for 
IRS 1, and in particular not for the low [Cii] intensity observed 
there. Solving the puzzle of the extremely low fine-structure in 
S 140 may be the puzzle to a general explanation of the fine- 
structure line deficit in galaxies with high infrared luminosities. 

At the cloud surface of S 140 we find a classical PDR, show¬ 
ing the expected stratification in the line tracers and producing 
strong [C ii] emission. The [C ii] line is narrower there, having an 
optical depth of f s; 3. The combination of the observed strat¬ 
ification and line intensities of [C ii] and [O i] can be explained 
by a plane-parallel PDR model with a density of 2.5 lO'* cm“^ 
if we assume that the surface is tilted by an inclination angle of 
80-85 degrees relative to the face-on orientation and illuminated 
by a a UV field of 60 xo, a value that is lower than obtained in 
the previous PDR model fits in the literature. 
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Figure A.l. Integrated map of [C ii] obtained in the GREAT ob¬ 
servations (contours) overlaid by colors of the [C ii] line intensity 
from the two PACS footprints. The top image shows the intensi¬ 
ties in the original PACS footprints. In the bottom version this is 
resampled to the positions of the GREAT map with an effective 
resolution of 17.2" beam. 


Appendix A: Comparison to Herschei data 

In the frame of the WADI key project dOssenkonf et alJl201 ih . 
S140 was observed with the HIEI and PACS instrument, also 
covering the [C ii] and [Oi] linefl. Part o f t he results were 
used e .g. in the papers of iDedes et al.l (l2010h and lKoumnia et all 
(l2015h . Due to the limited observing time in the key project only 
a small part of the source was covered. In the overlapping part 
we use the data here to check the consistency with our new, more 
extended observations with GREAT. 

A.1. PACS spectroscopy 

The [O i] and [C ii] lines were observed earlier in S 140 through 
the PACS instrument onboard Herschei in the frame of a spec¬ 
tral scan at two individual pointing, one centered at IRS 1 and 
the other one close to the cloud surface towards HD 211880. 
The PACS footprint is not fully-sampled so that it may miss part 

’ Herschei obsids: PACS: 1342222255 (IRS 1 [Oi] 63pm), 
1342222256 (IRS 1 68-208 pm), 1342222257 (Interface [Oi] 63 pm), 
1342222258 (Interface [Oi] 145 pm), 1342222259 (Interface [Cii] 
158 pm), HIFI: 134290781 ([Cii]) 



Figure A.2. Measured intensities in the different pixels from 
Fig. lA.lh as seen by GREAT and PACS. Green asterisks rep¬ 
resent the footprint around IRS 1, black plus signs the footprint 
at the outer interface. The brown line shows the identity, the blue 
line represent an identity with an offset of 30 K km s“’. 


of the flux between the pixels, but if we assume a source size of 
5.5" or above that missing flux should be small so that we can 
use the data to check the consistency with the GREAT observa¬ 
tions. The spectral resolution of PACS was insufficient to resolve 
the fine, so that we can only compare integrated intensities. 

Figure IATTI compares our new intensity map with the PACS 
data for the two footprints. We see that the spatial structure of 
the emission peaks around IRS 2 and the interface is consistent 
between PACS and GREAT data. The PACS data just miss the 
global emission peak, but they peak at the position closest to the 
true maximum and they also trace the structure of the extended 
emission at the outer cloud interface. When comparing absolute 
intensities, however, we find systematically somewhat lower in¬ 
tensities in the PACS data compared to GREAT data. 

To understand whether this is a calibration problem, we com¬ 
puted the pixel-to-pixel statistics in the resampled map. This is 
shown in Fig. IA.2I We find a uniform behavior between both 
footprints; the PACS intensity is offset from the identity by about 
30 K km s '. This means that we do not have a calibration dif¬ 
ference, but a difference in the absolute level. This can be nat¬ 
urally explained if there is a 30 K kms“* contamination of the 
OFF position used as the reference in the PACS observations. 
The GREAT observations use a well controlled OFF position 
6.5' south-west of the molecular cloud, while the chopped PACS 
observations involve a reference that falls 6' into the molecular 
cloud. Some [C ii] emission from that region is quite likely and 
would explain the constant absolute offset between the GREAT 
and PACS intensities. 

As the [O i] emission is much less extended, we expect no 
OFF contamination in the PACS data for this fine. Figure IA.3I 
compares tfie spatial distribution of the [O i] integrated fine emis¬ 
sion measured by PACS with our GREAT observations. The 
spatial distributions of the emission seen through GREAT and 
PACS agree relatively well, having the peak again around IRS2. 
But for the individual intensities we find some significant de¬ 
viations. Fig. IA.4I compares PACS and GREAT intensities from 
two different ways of resampling. One is shown in the lower 
panel of Fig. lAJl namely both maps are resampled to a common 
4" grid with 12" resolution (green asterisk and black plus sings 
in Fig. IA.4I) . Since the PACS map is not fully-sampled, we ex- 
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Figure A.3. Integrated map of [O i] obtained in the GREAT ob¬ 
servations (contours) overlaid colors of the [Oi] line intensity 
from the two PACS footprints. The top image shows the intensi¬ 
ties in the original PACS footprints with GREAT data convolved 
to 10" resolution. In the bottom version both maps are resampled 
to the positions of the GREAT map for the effective resolution 
of a circular 12" beam. 


tracted the flux based on the geometrical overlap. In the second 
method, we used the integrated intensity of the original 5x5 
PACS spaxels, and resampled the GREAT data to each spaxel 
position (blue diamonds and orange triangles in Eig. lA.3l In both 
ways the resulting trend is the same. There is a good match for 
the map around the interface and for the highest intensities close 
to IRS 2 at the northern boundary of the PACS held. However, 
across the PACS array centered at IRS 1, the intensities seem 
to by systematically too high. Off-contamination on the GREAT 
side can be excluded as it would show up as a constant offset 
like seen in the [C ii] case. Moreover, we used the same refer¬ 
ence position as for [C ii] and [O i] is rather less extended than 
[Cii] due to the higher critical density. On can see in Eig. IA.4I 
that the different ways of resampling the data do not make any 
significant difference. The PACS observer’s manual (v2.3) does 
not provide any direct explanation. 63 pm is not listed as a ghost 
wavelength from other strong emission lines and it is not in the 
leakage wavelength range. We can only speculate that this may 
be a stray-light problem where the strong [Oi] emission just at 
the edge of the PACS footprint is partially picked up in the ad¬ 
jacent spaxels. Unfortunately, there are no calibration data avail- 



FigureA.4. Measured intensities in the different pixels of 
Fig. lA.3l as seen by GREAT and PACS. Green asterisks represent 
the footprint around IRS 1 and black plus signs the footprint at 
the outer interface when both data sets are resampled to a com¬ 
mon 12" resolution (Fig. IA.3b '). Blue diamonds and orange tri¬ 
angles are obtained when we resample instead the GREAT data 
to the 25 pixels of the PACS array using the 9" PSF of PACS for 
the IRS 1 region and the interface respectively (Fig. lA.3h l. The 
brown line shows the identity. 

able to verify that situation as the calibration observations were 
not set up to just miss the target. Hence, we cannot provide a 
quantitative explanation for the mismatch. 

A.2. HtFt cut 

The [C It] line was also observed by HIFI (iDedes et alj2010h in a 
single strip going through IRS 1 and the surface also oriented in 
the same angle of 37 degrees east of north as our GREAT map. 
In Fig. lA.5l we compare individual HIFI spectra along that cut at 
a half-sampled spacing with the corresponding GREAT spectra 
obtained by averaging all data in a 6" circle around the HIFI 
positions. 

We And a generally very good match between HIFI and 
GREAT spectra, in particular the line shapes agree in all details, 
but at some positions the GREAT peak intensities are larger by 
a few up to 15 %. First of all this confirms the good relative cal¬ 
ibration and pointing accuracy of both instruments/telescopes. 
The higher intensity of GREAT at an offset of -50" and h- 12" 
could be due to the somewhat larger telescope beam of SOFIA 
leading to more intensity pickup from the surface at -50" and 
from the newly detected intensity peak north of IRS 1 at -1-12". 

Appendix B: The [Oi] 145 pm/63 pm ratio 

In spite of the possible stray-light issue of the PACS 
[Oi] map (see Appx. lA.ll l. we can also try to use the 
[Oi](145 /im)/[Oi](63 pm) line ratio to constrain the parame¬ 
ters of the PDRs in S 140, assuming that the ratio between both 
lines should be more robust than the simple integrated line in¬ 
tensity. 

Because the excitation energy of the 145 pm line of 
325 K lies above that of the 63 pm line (228 K) the 
[Oi](145 /rm)/[Oi](63 pm) measures temperature (or gas heat¬ 
ing through the U'V field) in the temperature range of about 
300 K. In many cases, however, this ratio may be dominated 
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Figure A.5. [Cii] spectra taken with HIFI (black) and GREAT (red) along a 37 degrees cut through IRS 1 and the cloud surface. 
The HIFI cut extends deeper into the cloud and only every second spectrum is plotted. The corresponding GREAT spectra were 
averaged in a 6" circle around the HIFI position. 
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of the [Oi] 145 pm line by a factor four, tow times lower than 
inferred from the [Oi](63 /tm)/[Cii] ratio in Sect. EJ 

At IRS 2, [Oi](145 /tm)/[Oi](63 pm) from PACS observa¬ 
tions is 0.3. For the PDR with a UV field of 5.6 x lO'^xo and a 
density between 1.4x 10^ and 10® cm~^ (Sect. l43T l. the predicted 
ratio is 0.03-0.045. Since the absorption by the cold gas at the 
back side of the PDR should affect only the 63 pm transition, a 
factor 7-10 of self-absorption is required to be consistent with 
the model predicted [Oi](145 jum)/[Oi](63 pm. 

At IRS 1, the observed [Oi](145 jum)/[Oi](63 pm) is 0.1, 
which suggests the PDR gas with a density of 10^-10^ cm~^ and 
a UV field of > 10®. But as discussed in Sect. 14.51 there is no 
consistent model for IRS 1. 


Figure B.l. Observed ratio of the [Oi] 145 pm to [Oi] 63 pm 
line integrated intensity seen by PACS (colors) overlaid on the 
contours of the [O i] 63 pm emission measured by GREAT. 


by the optical depth in the 63 pm line (iKaufman et al.l [T999t 
iMalhotrai et al.ll200T1) . A low [Oi](145 /im)/[Oi](63 pm) com¬ 
pared to the PDR model prediction is observed in many PDRs, 
and the foreground absorption, opacity effect of the [O i] 63 pm, 
and the cloud geometry are proposed as possible explana- 


_ etry _ 

tions (iLiseau et al.lll9^1caux et ai.lll999HGiannini et alj _ 

lOkada et al. 2003 ). lL'iseau et al.l (l2006t) moyided the theoretical 
estimates of this ratio for different physical conditions. In case 
of the optically thin emission, the ratio oyer 0.1 can be achieyed 
only when the gas temperature is > 500 K and the collision 
partner is H 2 , which is unlikely. For optically thick case, their 
calculation with the density of 3 x 10"^ cm“^ shows that the col¬ 
umn density of V(H)> 10^“* cm“^ would be required to make the 
[Oi](145 //m)/[Oi](63 /zm) ratio aboye 0.1. 


;xplana- 

Hbooot 


Figure ED shows the integrated intensity ratio of 
[Oi](145 jum)/[Oi](63 pm) from the two PACS footprints cen¬ 
tered at IRS 1 and at the int erface obseryed in the frame of the 
Herschel key project WADI dOssenkonf et al.ll201 ih . For the in¬ 
terface parameters discussed aboye, the face-on PDR model pro¬ 
duces a [Oi](145 jum)/[Oi](63 pm) ratio of 0.04, whereas the 
obseryed yalue is 0.17, i.e. about four times higher. If we inter¬ 
pret it through the inclination of the PDR in the same way as 
for the [C ii] line, the ratio asks for a geometrical amplification 









































